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1. Introduction 

Atmospheric neutrinos originate from the de- 
cay of unstable particles produced by cosmic ray 
interactions in the upper atmosphere. The use- 
ful spectrum runs from about 200 MeV to about 
1330 MeV. But flux extends to much higher ener- 
gies. The Super Kamioka group has studied these 
out to about 10 GeV. The higher energy sample 
is only about 21% of the size of the sample below 
1330 MeV. 

Atmospheric neutrinos approximate the ideal 
of a two distance neutrino oscillation experiment. 
Figure [l] shows the log of the flight length for 
neutrinos as a function of the cosine of the zenith 
angle. Events from below travel on the order of 
8000 km. Those coming from above travel on the 
order of 10 km. There is rather little solid angle 
at intermediate distances. Atmospheric neutrinos 
provide a good mixture of electron and muon neu- 
trinos since a significant number of the muons at 
these energies decay before reaching the ground. 
Simple arguments suggest that the ratio of muon 
to electron neutrinos should be about 2. 

To summarize, atmospheric neutrinos are a 
mixture of muon and electron neutrinos coming 
from hadronic and muon decays. The source 
is approximately at two distances that differ by 
three orders of magnitude. Most of the flux spans 
about a factor of 7 in energy. 

The atmospheric neutrino anomaly Q is the ob- 
servation of an apparent deficiency of muon type 
atmospheric neutrinos. The anomaly is believed 
to be evidence for neutrino oscillations But 
some of the details and assumptions used in the 
analysis cast some doubt on the oscillation hy- 
pothesis. The neutrino sample is being drawn 
from natural sources which adds a degree of un- 



certainty. 

Almost all analyses have assumed that the 
anomaly is restricted to the muon neutrinos and 
that the electron neutrino flux is a good normal- 



2. R(E y ) 

The ratio of observed to expected neutrino in- 
teractions is a measure of the possible reduction 
in the neutrino flux that is indicative of neu- 
trino oscillations. Due to a significant uncer- 
tainty in the overall flux normalization, as much 
as 25%, the quality of the attenuation statistic 
is improved by normalizing to something which 
is better known. The value of R = Hi 



which attempts to normalize the muon neutrino 
flux to the electron neutrino flux is thought to 
have significantly smaller systematic errors, on 
the order of 5%. 

An advantage that R has over other observ- 
ables is that it maximizes the use of the data. It 
is insensitive to details of the angular correlation 
between neutrino direction and reconstructed lep- 
ton direction from the interaction. A disadvan- 
tage is that one is averaging over two different 
baselines so that if only the long distance sample 
has oscillated significantly the result will be di- 
luted by the unoscillated sample from the shorter 
baseline. 

2.1. The Fit to R{E V ) 

Figure || illustrates the Super Kamioka mea- 
surement of R. The vertical scale is , 1 ^,"?„x which 
combines the observed data, R with a fit param- 
eter sin 2 (29). The horizontal axis of figure || is 
the logio of Am 2 /-E„ with Am 2 in eV 2 and E v 
measured in eV. The points at the center of the 
plot are taken from references [g] and |j| and are 
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Figure 1. Log of the flight distance in km as a function of zenith angle. Downward going is on the right. 
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Figure 2. si ^f^e) as a function 01 the l°gio of Am 2 /E u with Am 2 in eV 2 and E v measured in eV. (Add 
9 to get E v measured in GeV). The points plotted near the center of the figure are from references [5] 
and [6] using the Am 2 from reference [7]. The points plotted to the right are our fit to Am 2 using the 
data of [5] and [6]. 
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plotted with sin 2 (26) = 1 which is obtained from 
the Super Kamioka fit described in the next sec- 
tion. Since, for these point sin 2 (28) = 1 the 
value of 1 — R can be read off the graph. Ref- 
erence § quotes R = 0.61 ± 0.03 ± 0.05 for 
energies below 1330 MeV. Reference || quotes 
R = 0.66 ± 0.06 ± 0.08 for the higher energy sam- 
ple. 

The diamond points in figure are the same 
data but the values of sin 2 (28) and Am 2 are a fit 
to these data points. The fit yields sin 2 (2#) = 0.8 
and Am 2 = 0.24 eV 2 . The points near the center, 
the Super Kamioka fit, would agree with the data 
much better if the points could be lowered. But 
to do this would require sin 2 (28) — 1.5 which is 
unphysical. It is difficult to assign a confidence 
level to our fit. The points share many common 
systematic errors. One would need to quantify 
these correlated errors in order to achieve a le- 
gitimate confidence interval. Since we have plot- 
ted Am 2 /E^ the higher energy data points are 
to the left. The Am 2 mass scale is driven by the 
low value of R measured for the higher energy 
data points. A reasonable fit can be obtained for 
Am 2 > 0.1 eV 2 . More details of the fit can be 
found in ||. 

The curve plotted in figure || is an integral over 
the distances of figure [| for 1 — R. For single 
component neutrino oscillations the result is a 
function of Am 2 / E v and proportional to sin 2 (28) . 
The curve illustrates the structure expected from 
the two distance scales involved. At very small 
Am 2 all path lengths are small compared to the 
oscillation length and so the value of R is 1. At 
very large Am 2 the oscillation length is small 
compared to all path lengths in the problem so 
that 1 — R will be 0.5 x sin 2 (26>), For a range of 
intermediate values of Am 2 one path length is 
large compared to the oscillation length and the 
other is small so that 1 - R « 0.75 x sin 2 (28), 

2.2. Geomagnetic Effects 

Figure || illustrates the energy dependence of 
R for various values of Am 2 . The figure was 
obtained by numerically integrating the Honda 
97 muon neutrino fluxQ. For modest values of 
Am 2 the value of R is energy independent over 
the range indicated. But it can never drop much 



below 0.75. As Am 2 rises portions of R drop to 
below 0.75 and start to approach 0.5. Only at 
values of Am 2 of about 0.1 eV 2 is energy inde- 
pendence restored with R < 0.75. Figure ||| was 
generated for sin 2 (28) = 1, its maximum value. 
For smaller values of sin 2 (28) the allowed values 
of R are greater. This plot supports the fit of 
figure H in that the only curve that gives an en- 
ergy independent value for R(E V ) at a value of R 
below «0.75 is Am 2 = 0.1 eV 2 

The curve of figure was done assuming an 
isotropic flux. A check using the Honda 97 flux 
reproduces the results to about 1%, except at the 
lowest energy point, E=250 MeV, where the dif- 
ference between isotropy and the calculated flux 
is less than about 3% for Am 2 < 1 x 10~ 4 eV 2 . 

2.3. Physical Interpretation 

The physical significance of a low value of R, 
if interpreted in terms of only a muon deficiency 
due to muon neutrino oscillations integrated over 
energy and distance, is that both distance scales 
must participate. Since the short distance scale 
is about 3 orders of magnitude shorter than the 
large one the favored value for Am 2 must be com- 
parably larger. 

3. The Super Kamioka Fit 

3.1. Systematic Errors 

The Super Kamioka group has fit data to the 
neutrino oscillation hypothesis using: 

X 2 (sin 2 (28) , Am 2 , t) = V {NLta ~ ^c) 2 + y $ 

<Ji includes the statistical error for the observed 
number of events in the i'th bin. aj is the error 
on the j'th fitted "systematic" parameter, (jj is 
assumed to be known. e*is a vector of 8 systematic 
parameters. The sum is over 70 bins; 5 bins in 
solid angle, 7 bins in energy for both electron and 
muon like events. 

The number of Monte Carlo events is estimated 
from 

, T i £*data X ^ 

N 1 mc = 7 x 2^ w 

M C events 
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Figure 3. A plot of R{E V ) for sin 2 (26>) = 1 and various Am 2 . Am 2 runs from 3.2 x 1CT 4 eV 2 to 0.1 eV 2 
as marked in the legend. These correspond in general to the ordering of the curves from top to bottom. 
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Super Kamioka fit parameters 



5 



Where the C is the relative livetime of the data 
or Monte Carlo and 



'As 



cos(6)) 



w= (l + a)(El/E ) s (l + 



xf e ^(sm 2 {2d),Am 2 , (1 + \)L/E V ) 

{ (l-A/2) 
(l+ft/2) 
x < (1 - /? m /2) 
(l + /3 m /2)(l 
. (l + /3 m /2)(l 



p JVpc ' 
2 JV„ - 

2^ 



The meaning of the various a, S, rj, (3, A and 
p are given in table [l]. The choice of term from 
within the large braces depends on which portion 
of the data is being compared. 

There are a number of problems with this fit. 
In particular the expression (N data — N % MC ) 2 in- 
cludes a substantial systematic error in N l MC 
which is not included in the denominator a 2 . In- 
cluding this error would tend to lower the \ 2 - 

But the systematic error in N\ IC is common to 
many terms of the sum. In fact a definition of x 2 
including this common error would need to be of 
the matrix form: 



(Ndata-NMC) {O- ) '{N dat a-N M c) 



Where (Ndata — Nmc) is a vector of all of the 
differences and (c 2 ) -1 is the weight matrix, the 
inverse of the covariance matrix, that includes the 
correlated errors. 

Even with these definitions it would be hard 
to interpret the resulting \ 2 m terms of a confi- 
dence level since the a priori distributions of the 
systematic errors are not known. 

3.2. The Asymptotic Assumption 

The fit mentioned in the section above includes 
an unjustified assumption. It is assumed that the 
"systematic" parameters may be fit in an unbi- 
ased way. Many of these parameters are corre- 
lated with quantities of physical interest. For ex- 
ample the up to down flux ratio, ?y s is closely tied 
to any observed angular asymmetry. The muon to 
electron neutrino flux ratio, f3 s is essentially the 
denominator if the R value. The use of /3 S and 
P m effectively decouple the fit from the measured 
value of R. 



The fit assumes that the systematic parameters 
may be fit to some portion of the solid angle where 
the data sample is not anomalous. In essence the 
fit assumes that the upper hemisphere represents 
a good sample of unoscillated events. But as dis- 
cussed in section II the value of R makes this as- 
sumption rather risky. Any R interpreted solely 
in terms of muon neutrino oscillations must in- 
clude the upper hemisphere. On the other hand, 
R could be anomalously low because the electron 
sample is not as expected. 

The angular fit is relatively insensitive to a 
Am 2 in the range of 5 x 10~ 2 eV 2 to 10~ 4 eV 2 . 
This is because the data has been collected from 
primarily two baselines separated by 3 orders of 
magnitude. Am 2 's over a rather broad range will 
give comparable fits to the data. It is only the 
relatively small amount of data near the horizon 
that is sensitive to these intermediate ranges of 
Am 2 . In general it is easier to study such ranges 
by looking at the energy dependence after, per- 
haps, integrating over all solid angle, as we have 
done in the earlier section of this note. 

4. Time Variation of the Event Rate and 

R 

The event rate and the measured value of R 
for data collected during most of 1998 was sig- 
nificantly lower than that used in earlier Super 
Kamioka papers on atmospheric neutrinos [fL0|. In 
particular the overall event rate during the 321.6 
live days of this period dropped by 12±3% as 
compared to the event rate noted in the earlier 
414.4 day period. A careful comparison of these 
newer data indicates that the muon event rate did 
not change. The 1998 muon rate was 3±5% above 
the earlier rate. But both the electron sample 
and the multiprong event rate declined sharply. 
The electron sample was down by 18±5%. Table 
H summarizes the differences between the 1998 
sample the earlier sample. 

The change in the measured electron rate 
means that R must have changed. R rose from 
0.61±0.03 to 0.76±0.04. This rise in R brings 
it from a region which is not consistent with the 
hypothesis of neutrino oscillations over approxi- 
mately half of the solid angle to a value consistent 
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Cumulative [§ Phys. Lett.g 1998 Sample 



Single Ring 


3224 


1883 


1341 


e-like 


1607 


983 


624 


/i-like 


1617 


900 


717 


Multi-Ring 


1271 


784 


487 


Total 


4495 


2667 


1828 


n/e 


1.01±0.04 


0.92±0.04 


1.15±0.06 


R = {n/e) bs/(^/e) M c 


0.67±0.02 


0.61±0.03 


0.76±0.04 


Exposure (days) 


736 


414.4 


321.6 


Event Rate 


6.11±0.09 


6.44±0.12 


5.68±0.13 


Single Ring Rate 


4.38±0.08 


4.54±0.10 


4.17±0.11 


e rate 


2.18±0.05 


2.37±0.08 


1.94±0.08 


H rate 


2.20±0.05 


2.17±0.07 


2.23±0.08 


Multi-Ring Rate 


1.73±0.05 


1.89±0.07 


1.51±0.07 



Table 2 

Comparison of Super Kamioka Sub-GeV Samples 



with such a hypothesis. Unfortunately there is 
no good reason to doubt the credibility of cither 
data sample. Super Kamioka has combined these 
samples in their analyses. The oscillation analysis 
utilizes the electron events, which are not consis- 
tent through this period. 

R depends on the ratio of observations to sim- 
ulations. The simulations did not change in any 
significant way between these two periods ]To|] 

So there is clear evidence that the electron sam- 
ple does not provide a reliable normalization for 
the analysis of atmospheric neutrino data. The 
assumption that the ratio of atmospheric electron 
to muon neutrino flux can be calculated to an ac- 
curacy of 5% is in error, unless some portion of 
the observed signal is not of atmospheric origin 
or is not neutrinos. |Tc|| 

The original Super Kamioka paper on the at- 
mospheric neutrino flux(^] notes, "Using a de- 
tailed Monte Carlo simulation, the ratio ^/f>" ata 
was measured to be 0.61±0.03(s£a£)±0.05(sys), 
consistent with previous results from the 
Kaimiokande, 1MB and Soudan-2 experiments, 
and smaller than expected from theoretical mod- 
els of atmospheric neutrinos" . This is no longer 
true. The 1998 data is not consistent with previ- 
ous results of the Super Kamioka detector or any 
of the earlier experiments cited. An R value of 
0.76±0.04 seems to be unique to this period. 

The multi-GeV data sample has no evidence of 



a rate variation in 1998. 

Independent confirmation of the drop in event 
rate would help ally fears that it is due to instru- 
mental or systematic effects. The Super Kamioka 
group has maintained very good monitoring of the 
detector performance Jll| so such an explanation 
seems unlikely. 

5. Conclusions 

Several factors indicate that the atmospheric 
neutrino anomaly may be more complicated than 
generally thought. In such cases, it is important 
to seek corroboration of the results from other 
experiments. It is also important to be aware of 
simplifying assumptions that may have gone into 
the interpretation or the analysis of the data. One 
must be careful not to ignore information simply 
because it is in conflict with a preferred solution. 

It would appear that the electron neutrinos can 
no longer be used for relative normalization, at 
least until the variation is understood. In the 
early days the anomaly was considered to be ei- 
ther a deficiency of muon like events or an ex- 
cess of electron type events. In fact, the results 
looked a bit like a systematic error in the classifi- 
cation method used for electron and muon events. 
Since a random error would move events from the 
larger population, the muons, to the smaller pop- 
ulation, the electrons. The observed event rate 
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for all events was much closer to most flux models 
than might be expected based the errors quoted 
for the estimates. Studies using the independent 
muon decay signature corroborate the morphol- 
ogy based results for muon and electron classifi- 
cation. 

The "neutrino" events observed in under- 
ground detectors do not agree with expectations 
based on a simple source of atmospherically pro- 
duced hadron and muon decays. Uncertainties in 
these expectations does not provide much guid- 
ance as to what might be wrong. Simple one 
component muon neutrino oscillations are unable 
to provide a convincing explanation for the small 
value of R nor for recent apparent drop in the 
interaction rate. 
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